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I. INTRODUCTION 



The AdS/CFT correspondence [lH3|, or more generally gauge-gravity duality, demonstrates a deep and fascinating 
link between black hole physics and the plasma phase dynamics of certain (holographic) strongly coupled gauge 
theories. Over the last 10 years the gravitational side of the correspondence has also started to be used as a tractable 
theoretical model of strongly interacting non-Abelian media with properties similar to quark-gluon plasma studied 
first at RHIC and now also at the LHC (see P| for the most recent review of these developments). Early efforts in the 
applications of gauge-gravity duality methods to hot QCD matter were motivated by hydrodynamic simulations of the 
expanding fireball created in heavy ion collisions and focused on obtaining transport properties of holographic plasmas 
by analyzing low-lying quasinormal modes and linear response theory. These results provided concrete numerical 
predictions for the simplest transport coefficients of strongly coupled non-Abelian media with TV = 4 super Yang-Mills 
plasma as the primary examplj^ and have eventually led to the formulation of fluid-gravity duality [6] . Fluid-gravity 
duality is a correspondence which maps solutions of relativistic Navier-Stokes equations describing holographic liquids 
to long-wavelength distortions of black branes in higher dimensional geometry. The direct connection between the 
dynamics of black objects in higher dimensional spacetimes and solutions of nonlinear hydrodynamics provided an 
opportunity to understand black brane geometries and their features in terms of dual fluids, as well as to gain insights 
about hydrodynamics from the properties of Einstein's equations. These perspectives, as well as the possibility of 
applications, have generated significant interest in the nonlinear dynamics of black brane spacetimes. 

Dynamical black holes and their characterization has also been an important research theme in mathematical relativity 
for the last couple of decades (see 0-Q for useful reviews of this subject). The exact characterization of a dynamical 
black hole has proven to be a surprisingly thorny theoretical problem for general relativity. The standard textbook 
definition associates black hole interiors with regions of spacetime from which no signal can ever escape [lo| . Thus, 
finding the exact extent of such a region is necessarily a teleological procedure: properly defining "ever" and "escape" 
means that one must examine the ultimate fate of all signals from a point before ruling whether or not that point is 
part of the black hole. Identifying the event horizon boundary of a causal black hole is similarly teleological. Thus, 
even though an event horizon is a congruence of null geodesies obeying the same rules as any other congruence, its 
evolution can appear to be acausal. For example the area increase of an event horizon is not directly driven by 
infalling matter or energy; instead the actual effect of an influx through the event horizon is a decrease in its rate of 
expansion. 

These observations are not just mathematical curiosities. The non-local nature of the event horizon is acceptable as 
long as one treats it as a causal boundary removing the region containing a curvature singularity from the dynamics 
of the rest of spacetimjf] and does not associate any physical characteristics with it. However, this is not the only 
role of the event horizon - for the last four decades, one of the most celebrated results of black hole physics has 
been the link between the area of the event horizon and entropy. Already in the 1960s it was established that event 
horizons necessarily increase in area jlo| and this has usually been interpreted as being equivalent to the second 
law of thermodynamics. Thus, the apparently acausal expansion of event horizons would seem to imply a similarly 
acausal evolution of entropy. This leads to problems since the origin of black hole entropy needs to be sought within 
microscopic theories underlying gravitational interactions in the sense of the holographic principle [lli [l^ . For 



^ For an excellent review of these early developments see Q|. 



^ Hence guaranteeing consistency of low energy description of black holes in terms of classical gravity. 
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asymptotically flat or asymptotically de Sitter spacetimes such theories are not known (in principle one could imagine 
any of these being non-local [13]), but in the anti-de Sitter context these are local quantum field theories in a suitably 
understood large Nc limit. At the superficial level it is hard to judge whether acausality of the event horizon is a real 
problem in the first two cases, but in the AdS/CFT context it definitely is. 

The teleological nature of the event horizon is one of the main motivations for the ongoing research program to 
characterize black holes (quasi)locally, identifying their interiors from the presence of strong gravitational fields rather 
than on the basis of the causal structure of the entire spacetime. Quasilocal horizons go by such names as trapping 
[l^ . isolated or dynamical horizons (the last two are both reviewed in Q). In all cases though, these horizons can 
be thought of as generalizations of the classical apparent horizons [lo|. Recall that apparent horizons are associated 
with foliations of spacetimes. Areas of strong gravitational field arc identified with the region on each surface that 
is covered by trapped surfaces. The boundary of that region is an apparent horizon and it is the outermost surface 



for which the outgoing light front does not expand in areslj. With a slight abuse of terminology the union of such 
surfaces over all time slices is often also referred to as an apparent horizon and it is in this sense that it will be used 
it here. 

In the context of gauge-gravity duality characteristics of black holes with planar horizons in higher dimensional 
spacetimes are at the same time the quantities describing dual finite energy density or finite charge density states 
of local quantum field theories. In static situations, where the acausal nature of the event horizon plays no role, 
the entropy defined by the event horizon was identified with the thermodynamic entropy of a dual holographic field 
theory. Such entropy satisfies a very strong constraint: the first law of thermodynamics linking IR quantities (i.e. 
temperature and entropy density) with UV quantities (the energy density) . Since the latter are well defined in the dual 
quantum field theories (energy density is the expectation value of one of the components of the energy-momentum 
tensor in thermal equilibrium), there are no controversies with associating thermodynamic entropy with the event 
horizon in time-independent situations. However in static situations (at least in the context of Kerr-Newman black 
holes/branes) the event horizon and one of apparent horizons coincide, so that by associating the entropy with the 
event horizon one actually associates it at the same time with an apparent horizon. The latter identification actually 
turns out to be more robust, as the example of conformal soliton flow [3l suggests [l^]- 

Beyond equilibrium three problems arise. The first follows from the aforementioned nonlocality of the event horizon, 
the second comes from foliation dependence of apparent horizons, whereas the third one is related to the various 
ways in which one can associate points on any of horizons with points on the boundary (this is referred to as the 
bulk-boundary map [l^I)- The last of these is necessary to localize the entropy production in the dual field theory. To 
illustrate that the first issue is a serious problem quite disconnected from any ambiguities of the bulk-boundary map, 
one can consider the example of gravitational dynamics with a sharp distinction between a dual equilibrium regime 
without entropy production and a dynamical transition with dissipation. The relevant backgrounds, much in the 
spirit of the Vaidya solution, appeared in the context of the thermalization problem of strongly coupled non-Abelian 
media and describe gravitational processes in which the dual quantum field theory undergoes a transition between 
two equilibrium states in a finite time interval [isl . [lol . In such a situation the event horizon evolves past the bulk 
lightcones spanned by the transition region on the boundary. The latter patch of bulk spacetime is dual to the boundary 
region where the holographic field theory is in equilibrium, so that its thermodynamic entropy stays constant. This 



^ It is also required that ingoing light front shrinks in area and that inside this surface there are other surfaces, such that both ingoing 
and outgoing light fronts emitted from them shrink in area. For a more precise definition see Section IIIII or the review articles 0, 0] • 
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result strongly suggests that the causal boundary of a black hole is not the relevant entropy carrier regardless of any 
ambiguities of the bulk-boundary map In contrast with the area of the event horizon, the entropy defined by 

the unique apparent horizon respecting symmetries of 1-dimensional boundary dynamics considered in [isl . [3| was 
constant before the transition process and eventually in the far future agreed with the one given by the event horizon. 

In less symmetric situations, apart from the choice of bulk-boundary map, the foliation dependence of apparent 
horizons becomes a significant issue - different foliations of spacetime lead to different apparent horizons. The 
most trivial and at the same time pessimistic possibility is that the notion of local entropy does not extend beyond 
equilibrium situations and foliation dependence, as well as the freedom of bulk-boundary mapping, signal exactly 
this. It might also be that on the dual field theory side there are many relevant local notions of entropy and different 
apparent horizons correspond to such different notions. Yet another possibility is that the field theory notion of 
entropy does not suffer from ambiguities of kinds introduced by foliation dependence of apparent horizons, which 
might be used as a guiding principle for finding preferred apparent horizon. 

Resolving these issues in the general case is very difficult if not impossible, so the only hope is to proceed example 
by example. In global equilibrium the foliation dependence essentially vanishes and the event and apparent horizons 
coincide. In the near-equilibrium regime one expects the horizons to be "close" (in a sense of [2^ or 21|) so that a 
dynamical apparent horizon behaves almost like an isolated horizon. If one takes the near-equilibrium regime as point 
of departure for further studies, one is immediately led to considering apparent horizons in the geometry of fluid- 
gravity duality. This background captures the hydrodynamic regime on the field theory side starting from a locally 
boosted and dilated black brane supplemented with gradient corrections 0] and from this perspective hydrodynamics 
can be regarded as the simplest (because of its universality) type of collective dynamics that quantum field theory 
can undergo. 

The generalization of entropy to hydrodynamics is provided by the notion of an entropy current. Such a current is 
constructed phenomenologically in the gradient expansion by requiring that in equilibrium it reproduces thermody- 
namic entropy and that its divergence evaluated on solutions of the equations of hydrodynamics is non-negative. A 
detailed analysis of the consequences of this generalized second law of thermodynamics on the form of the entropy 
current in [22] showed that even up to second order in gradients there is an ambiguity inherent in such a definitioro. 
From the point of view of fiuid-gravity duality it was very natural to ask what is the gravity interpretation of the 
coefficients appearing in the boundary hydrodynamic entropy current and what might be the bulk counterpart of the 
ambiguity in its definition. In the pioneering work 17] a candidate entropy current was obtained by mapping the area 
theorem on the event horizon onto the boundary along ingoing null geodesies. In general, there are infinitely many 
directions in which such geodesies can propagate from the boundary, but hydrodynamic covariance requires that such 
geodesies close to the boundary move in the direction specified (at leading order of the gradient expansion) by the 
local fluid velocity [l^]. Ambiguities appearing in such procedure appear at third and higher orders in the gradient 



expansion and were irrelevant in the second order construction of |17l . |22| . This causal bulk-boundary map can be 
supplemented with suitably understood boundary diffeomorphisms and the latter turn out to capture precisely the 
ambiguity discussed in j^j. 

rn 

Furthermore, in a recent paper ^23] it was shown that for a fixed bulk-boundary map the same freedom in entropy 
current can be understood as coming from different bulk hypersurfaces with a fixed foliation satisfying a generalized 



Part of the ambiguity is trivial and comes from a term whose divergence vanishes. 
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area theorem and asymptoting to the event horizon. Such surfaces were dubbed "generahzed horizons" with the 
event horizon and the (asymptoting to it) apparent horizon being just two particular instances of the more general 
notion. From the perspective of the phenomenological definition of the hydrodynamic entropy current none of these 
hypersurfaces and none of the available bulk-boundary maps is favored over any other. However, causality of the 
boundary field theory seems to favor the entropy current dual to the apparent horizon - providing that it is free of the 
ambiguities related to foliation dep endence and that the bulk-boundary map in use is causal. The task of this paper 
is to elaborate on the proposal 2J| by presenting a derivation of the apparent horizon in the geometry of fluid-gravity 
duality, its features, as well as discussing the properties of the dual entropy current. 

The organization of the paper is the following. Section presents in a self-contained fashion the geometry dual to 
conformal fluid dynamics in arbitrary dimensions obtained in [25! . Section IIIII is the main part of the paper and 
provides the detailed calculation of the relevant apparent horizon in the case of conformal fluid-gravity duality up to 
second order in gradients. Section HVl focuses on the hydrodynamic side of fluid-gravity duality and analyzes the dual 
entropy current using the technology introduced in [2J|. The general discussions of the results and possible future 
directions of research are provided in Section |Vl Appendix |X] provides some details on the Weyl-covariant derivative 
and Weyl-covariant hydrodynamic tensors, whereas Appendix [B] illustrates the methods developed in Section HTll bv 
describing the construction of an apparent horizon in the Vaidya spacetime. Readers interested mostly in the general- 
relativistic aspects of these considerations can skip Section IIVI and regard the paper as an example of a perturbative 
calculation of an apparent horizon in a geometry governed by Einstein's equations with negative cosmological constant. 



II. THE GEOMETRY OF FLUID-GRAVITY DUALITY 



The geometry of fluid-gravity duality in arbitrary dimensions [25| is a solution to Einstein gravity with negative 
cosmological constant 



N 



d'^+'^xV^{R + d{d-l)}, (1) 



where Gn is the (d -I- l)-dimensional Newton's constant and the AdS radius is set to 1. The action ([T|) arises in the 
context of string theory (for d = 2,3,4 and 6, see [3] for details) and describes a sector of decoupled dynamics of 
the one-point function of the energy-momentum tensor operator in planar strongly coupled holographic conformal 
field theories 25|. The equations of motion derived from ([T]) support a d-parameter family of exact, static black hole 



solutions with planar horizons obtained by boosting and dilating the AdS-Schwarzschild black brane solution 

ds^ = 2'«^da;''dr - ( 1 ^ ) Ufj.u^dx^dx" + [g^i, + u^-u^) dx^dx" . (2) 

Here g^j/ denotes components of the flat Minkowski metric on the conformal boundary of the asymptotically AdS 
spacetime ([2]). The boost parameter is a d-component velocity in the directions, normalized so that u^u^^ = — 1 
in the sense of the boundary metric (7^,^. The lines of constant in 1^ are ingoing null geodesic, for large r 
propagating in the direction set by u^, and the radial coordinate r parametrizes them in an affine way [17]. The 
geometry ([2]) may be regarded as a stack of constant-r d-dimensional planes, starting from the boundary at r = cxi 
(which is d-dimensional Minkowski spacetime) right down to the curvature singularity at r = 0. The latter is shielded 
by the event horizon at r = l/b, which is at the same time an (isolated) apparent horizon. The dilation parameter b 
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appearing in ([2]) is related to the Hawking temperature T of the event horizon by 

"-^^ (=» 

Unhke black holes in asymptotically flat spacetime, the metric © supports perturbations varying much slower within 
the transverse planes than within the radial direction. The parameter controlling the scale of variations in the radial 
direction is 6. 

If h, and g^v are allowed to vary slowly compared to the scale set by 6, the metric ([2]) should be an approximate 
solution of nonlinear Einstein's equations with corrections organized in an expansion in the number of gradients in the 
directions. Direct calculations Q have shown that proceeding in this way is a systematic way of solving Einstein's 



equations, provided that the dual energy-momentum tensor [271 128[ depending on h and is conserved. For the 



metric ([2]) the dual energy-momentum tensor is that of a relativistic perfect fluid with a conformal equation of state 

T^"' = eu^u" + pP^"" , (4) 

where 

Piiv = gfiu + u^Ui, (5) 
is the projector operator onto the space transverse to , 

. = (.-l)p=(d-l)^5- (6) 
and g^i^ is some weakly curved metric in which the fluid lives. 

Both the metric ([2]) and the energy-momentum tensor (|4|) receive gradient corrections. The separation of scales 
mentioned earlier implies that corrections to the metric ([2|) will be tensorial quantities made of x'^-derivatives of 6, 
and g^u with scalar functions of r and h as coefficients. The relevance of these terms is suppressed by the number 
of gradients and for practical reasons the expansion is terminated at the 2-derivative level. The tensorial quantities 
in question are scalars S, transverse (u^F'' — 0) vectors V^^ and transverse {u'-''T^i, = 0) traceless symmetric rank 2 
tensors T^Jfl- A priori one should consider all possible terms, as was done originally in This task can however be 
greatly simplified by utilizing the underlying conformal symmetry and seeking Weyl-invariant solutions of Einstein's 
equations, i.e. solutions invariant under simultaneous rescalings of 

^ e-^-^g^^, w^^e'^M'', b ^ e'H and r ^ e'^r (7) 

where depends on the coordinates |25|. The leading order metric ^ is Weyl-invariant, but due to the presence 
of dr it does not retain its form at higher orders. It can however be written in a manifestly Weyl-invariant form upon 
introducing a vector field Av defined by [2^ 

Au = u^VxUi, — -Uu . (8) 

a — 1 

This quantity is of order one in the gradient expansion and transforms as a connection under Weyl-transformations 

A^^A^ + d^(t> . (9) 



It needs to be borne in mind, that all those quantities are also required to be independent when evaluated on lower order solutions of 
hydrodynamics . 
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The Weyl-invariant form of the metric ([2|) reads 

ds^ = 2M^da;^ (dr - Ar^dx'') - ( 1 ^ ) u^u^dx^'dx'' + (g^,^ + u^,u^) da;''da;''. (10) 

V {rh) J 

This metric is a leading order approximation to a spacetime whose metric is of the form 

ds^ = Gabdx'^dx'' = -2u^dx^(dr + V„dx") + g^^dx^dx" (11) 

with the condition u'^Gp.i' = completely fixing the gauge freedom. The subleading corrections to (flUl) need to be 
Weyl-invariant and the simplest way to construct them is by summing individual Weyl-invariant contributions order 
by order in the gradient expansion. A single Weyl-invariant contribution to (llip can be represented as a scalar 
function of the Weyl-invariant combination rb multiplying a Weyl-covariant (i.e. transforming homogeneously under 
Weyl transformations of b, and g^y, see Appendix \K\ tensor of a given weight w supplemented with a factor of h'" . 

A powerful tool in generating Weyl-covariant gradient terms is the Weyl-covariant derivative I?^, which uses the 
connection ([5]) to compensate for derivatives of the Weyl factor coming from derivatives of Weyl-covariant tensors. 
It has the property that a Weyl-covariant derivative of a Weyl-covariant expression is itself Weyl-covariant with the 
same weight (see Appendix 1X1 or the original publications [itI, 0, 29 1 for details). 



At first order in gradients there is only a single Weyl-covariant term available, which is the shear tensor of the fluid 



a^v. It reads 



^i^y = o Pp"<> + T^vUf,) (12) 



and transforms with Weyl-weight 3. At second order in gradients, there are in total 10 Weyl-covariant terms: 3 
scalars, 2 transverse vectors and 5 transverse traceless symmetric rank 2 tensors. For convenience these objects can 
be defined with appropriate powers of b to render them Weyl-invariant. The scalar contributions read 

Sl^b^af,^(J^"' , S2^b^io^^oj'"' and S3 = b'^TZ, (13) 

where uj is the vorticity of the flow and TZ is the Weyl-covariant curvature tensor and curvature scalar (see Appendix 
\X\ for details) . The Weyl-invariant transverse vectors are 

Vi^ = bP^,Vpa''P and V2 = bPf^^VpCu'^P (14) 

Finally, the Weyl-invariant tensors read 

Ti ~ u'"D pa , T2fiy ~ Cpai'pu°'u'^ , Ts^^ = w^o-p^ + w^Cp^, 
T4f_,u ^ (jPapy - -j^-jPfj^^aafia"''^ and Tsp^ = w^Wp^ -I- -^-^Pp^Wa^w"'^, (15) 

where Cp^avp is a Weyl-covariantized curvature tensor (consult Appendix |X] for its detailed form). 

The metric ([TT|) up to second order in gradients can be expressed in terms of (fT^ . (IT51) . and (IT5|) and takes the 
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form 



Vf, = r^Bu^ + rAf, 



2r 



(6r) 



5o 



1 



2(d- 1) {brY 



K2S1) 



^ Pn.b-^Si)-^K,SiP^ 



d-1 



(16) 



(17) 



where 



B 



2{br)' 



r(l-(H')- 



The quantities P, i?2, -f'^i, ^2, L are functions of 6r introduced in 25| and read 



(18) 



(19) 



roo d-2 _ I 



(20) 



H2{br) 



br 



a^' - 1) Ji 



y^-^dy [l + {d- l)yF{y) + 2y^F'{y)\ 



dy, 



(21) 



Ki{br) 



OG 1^ pOO 



br 



y^F'iyf 



(22) 



i^2(fer) ^ e(e - l)P'(e) - 2(d - 



(23) 



L{bT) = / e'^-^e / dy 



(24) 



The metric given above is a solution of Einstein equations with negative cosmological constant up to second order 
in gradients, provided that h and u^^ satisfy the equations of dual hydrodynamics, i.e. the equations of covariant 
conservation of the energy-momentum tensor obtained from (|16p by holographic renormalization 



(25) 



where 
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The geometric picture emerging is that of spacetime locally approximated by tubes of uniform black branes spanned 
along ingoing null geodesies given by lines of constant . The dilation and boost parameters b and u^, as well as the 
boundary metric g^u vary from tube to tube, but, as anticipated, the scales of these variations are small compared 
to variations of the bulk metric along the radial null direction [l^]- Due to this tubewise approximation, the leading 
order geometry of fluid-gravity duality inherits the causal structure of static black brane, i.e. the event horizon located 
at r = 1/h, now with h depending on a;^ [17]. It is to be expected, and is confirmed by direct calculation further in 
the text, that the event horizon of ([2]) with slowly varying h, u'^ and g^^ is at the same time an apparent horizon. 
Such an apparent horizon is called isolated and does not lead to entropy production. The isolated apparent horizon 
at r = 1/& is expected to become dynamical once corrections to ([2]) are included and its position will also be modified. 
The dynamics of this almost isolated apparent horizon can be described in a gradient expansion much in the spirit of 
the framework of slowly evolving horizons 30 -3^. 



III. LOCATING THE APPARENT HORIZON IN THE GEOMETRY OF FLUID-GRAVITY DUALITY 



This section is devoted to identifying an apparent horizon for the spacetimes defined by the metric (|lip . This search 
will be based on two criteria: 1) from the isolated horizon contained in the unperturbed geometry ([2]) it is natural to 
expect the apparent horizon to be a perturbation of the hypersurface r — l/b[x) and 2) to ensure compatibility with 
the dual conformal fluid solution those perturbations are required to be manifestly Weyl-invariant. 



A. Preliminaries 



Apparent horizons are defined in terms of trapped and marginally trapped surfaces. In both cases the term "surface" 
means a codimension-two hypersurface VL embedded in a larger spacetime. The normal space to such a surface is 
spanned at any point by a pair of null vectors t and n. The following considerations apply to spacetimes where it 
makes sense to specify that both of these are future-oriented and respectively outwards and inwards pointing. It is 
convenient to cross- normalize them so that ^ ■ n — ~1 (this leaves a degree of scaling freedom) . 

The induced metric on 17 can be written as 

qab = gab + Lnb + , (27) 
while the outward and inward null expansions of fl are 

Oii) = r'Va^ = Ce log ^ and 0(„) = q^'^aUb = £„ log (28) 
or, more generally, for an arbitrary normal vector X = A£ + Bn 

Oix) = A9(^,) + B9^,,) . (29) 



Now Q is said to be outer trapped if 9^£^ < 0, trapped if < and < and untrapped if = and < 0. It 
is outer marginally trapped if = and marginally trapped if 6'(£) = and 6'(„) < 0. Trapped surfaces are indicative 
of black hole regions, with well-known theorems linking them to both singularities and the existence of event horizons 



10| . As recalled in the introduction they are also used to define apparent horizons 10| . Given a foliation of spacetime 
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into spacelike hypersurfaces Et ("instants" of time) one can define the total trapped region on each St as the union 
of all the outer trapped surfaces. Then (up to some technicalities which will be ignored here) the boundary of each 
of those regions fit is outer marginally trapped and known as the apparent horizon. A common abuse of terminology 
(adopted in the following) also uses the term apparent horizon to refer to the hypersurface defined by the evolving Qt 
(that is the union of the fit). 

In practical calculations, this definition of an apparent horizon is not very usable and instead one just searches directly 
for hypersurfaces foliated by outer marginally trapped surfaces. This is common practice in numerical relativity (see, 
for example, 0] and references therein). More generally, the teleological nature of classical black holes and their 
event horizons has lead many to search for a (quasi)local and properly causal replacement. Horizons foliated by 
(outer) marginally trapped surfaces which (hopefully) bound regions of trapped surfaces are the most obvious and 
mathematically tractable candidates. 

For example, the boundaries of stationary black holes (or branes) are taken to be weakly isolated horizons: 
codimension-one hypersurfaces that are foliated by outer marginally trapped surfaces or isolated horizons if their 



33, 



34| or review articles such as [7h9|). These 



extrinsic geometry is also invariant (see for example the discussions in 
are closely related (though more general than) Killing horizons and under many circumstances do a good job of 
characterizing a stationary black hole boundary without reference to causal structure or infinities. This is particularly 
so if one adds extra conditions to ensure that there are fully trapped surfaces "just inside" the horizon. For Hayward's 
[l^ future outer trapping horizons (FOTHs) one assumes that 

6'(„) < and £„6i(<>) < . (30) 

That is, the inward expansion is negative and under a small inwards deformation the outward expansion also becomes 
negative. The black branes considered in this paper are examples of FOTHs. 

For the classical definition, it is clear that time-evolved apparent horizons are foliation dependent: different foliations 
will sample a different set of trapped surfaces and so give rise to a different "time-evolved" horizon. Alternatively, 



focusing on the time-evolved horizon itself, it can be shown (see, for example [31|,|35[) that a hypersurface foliated 



by outer marginally trapped surfaces is not rigid and may be deformed while maintai ning its properties. The non- 



uniqueness of apparent horizons has been explicitly demonstrated in several papers |36l . |37 1 



B. Finding the horizon: strategy 



Problems with uniqueness are somewhat alleviated in the present calculation by the xxrequirement that perturbations 
of the horizon be manifestly Weyl covariant. Then, the time-evolved apparent horizon A should be specified as the 
level set of a scalar function 



S{r,x) = b{x)r - g{x) , 
where g{x) is a Weyl- invariant scalar defined by 

g{x) = gi{x)+g2{x) + . . 



(31) 



(32) 
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where gk denotes a linear combination of all Weyl-invariant scalars at order k in the gradient expansion. There are 
no Weyl-invariant scalars at order 1, and 3 at order 2, so one expects to find 

31 (x) =0 

92{x) = hiSi{x) + h2S2{x) + h-iSsix) , (33) 

where the Si are the 3 independent Weyl-invariant scalars (jl3p and the constants hi will be determined by solving 
= and the conditions (|30|) . Once this is done, the expression for the position of the apparent horizon will take 
the form 

r = rnix) = 7 (1 + hiSi{x) + h2S2{x) + h^S^{x)) . (34) 


This is a strong constraint, but a reasonable one to impose in a perturbative regime where physical considerations 
suggest that the horizon should be given by a Weyl covariant structure. Testing these surfaces as potential horizons 
means that one must consider their possible foliations and find out whether any of them satisfy = and the 
conditions ([30]) . Again however, one can lean on the Weyl covariance to simplify the calculation. Specifically, the 
outer marginally trapped surfaces of A will have their own (in A) normal v. This vector is required to be expressible 
as a sum of Weyl invariant terms and further that it be surface forming 

1; A d^; = . (35) 

Though this only really needs to apply on the horizon itself, it turns out to be computationally much easier to check 
this condition for v specified not only on the putative horizon but also in some neighbourhood. Thus, in practice one 
should look for Weyl-covariant one- form fields that are surface forming in some neighbourhood of A. 

Thus the search domain will not be arbitrarily large, but rather be restricted to potential horizons and foliations that 
are essentially Weyl-covariant perturbations of the unperturbed boosted black brane solution ([2]). Marginally outer 
trapped surfaces are to be sought among intersections of these classes. It will be shown in the following that for 
the geometry of fiuid-gravity duality, up to second order in the gradient expansion the conditions =0 and (|35p 
determine v (as well as the hi in ()34p ) uniquely. 

C. Finding the horizon: hypersurfaces and intersections 

The program outlined above can be implemented as follows. The normal covector to a surface of the form ([3T|) is 

m = dS, (36) 

which up to second order in the gradient expansion is 

m = rdf,b dx^" + b dr . (37) 

The function g does not contribute above, since the leading term involves 9^(72, which is of third order in gradients. 
It is convenient to write the normal in terms of the Weyl-covariant derivative, which acting on b (Weyl weight —1) is 



'D^b = df,b - A^b 



(38) 
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One then ha.^ 

m = r{Vf,b + bAf,)dx>' + bdr. (39) 



Raising the index using the metric (jlip one gets (using the formula for the inverse given in 25|), up to second order 



+ 1 (d^ + Iwr^^' + (-d(d^^" - + Jdh^JhY^^'^^' + (Z3^ (d3T)^3) . (41) 

Next one must consider potential foliations of A. As noted earlier, the foliation of the apparent horizon can be 
specified by a vector field v which is tangent to A and but otherwise normal to the leaves. Given a parametrization 
of the horizon so that S{r{y), x{y)) = const, this means that 

dS dr d^dx^_r. 

In terms of this coordinate system the tangent vectors are, of course, d/dy" which push-forward into the full space 



time as 



and applying (|42|) one finds 



^ _dx^ i _d_ _ fdS_Y'^ ^ 



(44) 



By construction these vectors all satisfy v ■ m — (as they should). 



It is very convenient to choose the coordinates on the horizon — x'^. This should be reasonable as long as the 
horizon does not "fold over"; this should be the case in this perturbative, gradient expansion limit. This choice also 
has the advantage of making the bulk-boundary map trivial (as discussed in Section llVp . Then the tangent vectors 
to A can be written as 

and a general vector field tangent to the horizon is given by 



- (46) 



In terms of the coordinate basis in the bulk one then has 



Recall now that the geometry found in '25] satisfies Einstein equations provided that the equations of hydrodynamics are satisfied by 
the quantities b, u^^ (in terms of which also yl^ is expressed). These equations imply [2^ that Dfe is of the second order in gradients. 
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Requiring that the vector v be Weyl covariant fixes (up to second order) 



v^' = bu^ + {Vi^b'^ci + V2''b^C2 + u^'iSibei + ^2^62 + Sgfoeg)} (48) 

2 

d(d~ 1) ' 



v"^ ^ -rA^buf" + Sir— — - , (49) 



where ci, C2, ei, 62 and 63 are some constants. It is computationally convenient to normalize v so that 

+ i;2 = , (50) 

in which case the coefficients of the longitudinal terms in (|48|) vanish 

ei = 62 = 63 = 0. (51) 

The remaining coefficients (ci, C2) appearing in v are also not arbitrary. As discussed earlier, to ensure that the vector 
V defines a foliation one has to impose the Frobenius condition (I35|) . There are two types of terms, which turn out to 
be given by 

V[f_cd^Vp] = 

up to terms of higher order in the gradient expansion. This determines the coefficients ci and C2 

2 



Cl 



d{d-2) 



C2^^. (53) 
This way one finds that the foliation vector v is completely determined once A is fixed 

Since the Frobenius condition was imposed for the full spacetime (rather than just on the horizon), this vector v 
actually gives rise to foliation of the full spacetime, at least in a neighborhood of the horizon. 

It is interesting that one gets a unique result. It seems plausible that this will also be the case at higher orders in the 
gradient expansion. To see this, note that at a given order fc, v is entirely specified in terms of its v'^ components, and 
its component does not depend on the k-th order contribution to v^. In complete analogy with the second order, 
at order k will be a linear combination of all available transverse and longitudinal vectors. The hypersurface A is 
specified as the level set of a scalar function S{r, x) (see (I5T]) ). which at order k contains all the available hydrodynamic 
scalars of order k. The vector m normal to A is defined by to = dS, so the construction does not introduce any 
further coefficients to be determined. Now consider the normalization condition (|50p and expand the f ^ contribution 

gt^uvf^V - 2u^V^v''v'' - 2u^v''v'' + TO^ = 0. (55) 



In order to evaluate the fc-th order contribution to ([551) from v^^ it is sufficient to take the zeroth order metric. Note 
however that since at leading order is proportional to and CJ^^ is transverse, the first term on the left hand side 
of (|55p vanishes for all r. Since does not receive corrections from the k-th order w^, the only term which depends 
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on this is actually u^V^v^v^ . But since is also proportional to at leading order, the whole left hand side of ([55)1 
at order k depends only on the longitudinal contributions to at this order. If so, formula (|55p fixes them uniquely 
and does not constrain the transverse contributions. What is left at order k are scalar contributions to S'(r, x) and 
transverse contributions to . But at a given order, transverse and longitudinal quantities are independent, so the 
scalar condition 9i at order k fixes all the contributions to S{r, x). The transverse components of w^, relevant for the 
foliation of A, are likely to be fixed by the Frobenius condition psp in analogy with what happens at second order. 
It can be checked that the contributions in question will appear in the Frobenius conditions, but it seems difficult 
to show that by choosing the transverse parts appropriately one can satisfy Frobenius conditions at any order. One 
argument that this is indeed the case is that such a condition must be satisfied for on the event horizon at arbitrary 
order and in this case is fixed and given by m^. It would certainly be interesting to make these statements more 
precise. 



D. Horizons 



Dynamical quasilocal horizons are spacelike and so m should be timelike and v spacelike. Without loss of generality 
one can assume that m is future pointing and v is outward pointing. Then the null normals to the surfaces of constant 
S and V are 



v = (i-Cn 
m = £ + Cn , 



(56) 



where the scalar C is called the evolution parameter |30|, |3l| . In this case 



(57) 



The sign of the evolution parameter indicates whether A is spacelike or timelike (or null if C = 0). The signs of the 
coefficients in (j56p have been chosen to ensure that both £ and n are future-pointing, and £ is outward-pointing while 
n is inward-pointing. 

The null normals are then 



+ (^^'((^ (5^ + '-^^ + ^^''(-2 + (Z^ w + 2^ (ip 



+ ci 



1 1 1 



'd-2 4 



1 1 1 



d-2 2{d-l){br 



2{d-2) (d- 1) 



(58) 



(59) 



and 



2Sr2 \ ' 'drb {d-2) (r6)2 (rb) 
1 

~b' 



2^L + c,) + V2''{- ^ ^ 



{d-2) {rbf 



C2, 



(60) 
(61) 
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As a check on the results obtained so far one can determine the location of the event horizon and compare with [25| . 
Using the explicit form of the evolution parameter C (obtained from ((57)) ) 

it is easy to calculate where the event horizon is located. On general grounds one expects a result of the form p4l) . 
Solving C{rEH) = one finds ([M)) with 

u(EH) ^ 2{d^+d-A) 1 

1 d^{d~l){d^2) d{d-l) ' 



h. 



d + 2 



(EH) 

2 2d{d - 2) 

(EH) f_ 

^ d{d~l){d-2) ' 



hr"'=- ,,, ,,,, (63) 



which matches the results of 



25 



Note that because the event horizon is null it must be the case that v is proportional to m. Using (|63p and the explicit 
form of m and w, it may be checked directly that this is indeed the case. 

To determine the position of the apparent horizon one needs to calculate the null expansions from the forms ()28p 

0(1) = g"'Va4 and 0(„) = r'Vanh , (64) 

where the metric induced on the foliation slices is calculated from (j27p . Using the results of the previous section one 
finds (up to second order 



^(^) = ^'-'^ {'"'^^Fr {'^^-Jdh) + 2(^(6^^^ - ^^^^^^^1(1)) + 



(Z^ + \wr^ + ^ (^-1)^2) ^^ ) ) 



d-1 1/1 



br br \(br)^ 



S2+rK{Si] . (66) 



Note that the results are manifestly Weyl-invariant. In particular, there is no correction at first order (as required by 
Weyl invariance). With these results in hand, it is straightforward to determine the location of the apparent horizon 
by solving 0(£)(rA_f/) = 0. One again finds (p4| with 

^(AH) d + 2 



2 2d(d-2)^ 

(AH) ^ l_ 

d{d- l)(d-2) 



/^r' = - ... (67) 



Only hi differs from the result for the event horizon [25 1 



This computation is fairly lengthy. 
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The expression 



explicitly shows that the apparent horizon hes within (or coincides with) the event horizon in the sense that an ingoing 
radial null geodesic will cross first the event horizon and only then the apparent horizon, since r is an afRne parameter 
on such geodesies. It is also easy to check that the apparent horizon is spacelike or null 

as required. 



IV. THE HYDRODYNAMIC ENTROPY CURRENT DEFINED BY THE APPARENT HORIZON 

In hydrodynamics the entropy current is a phenomenological notion constructed order-by-order in the gradient expan- 
sion starting from the term describing the flow of thermodynamic entropy. Subleading contributions are given as a 
sum of all available hydrodynamic vectors (not necessarily transverse) chosen in such a way that the divergence of the 
current is non-negative when evaluated on the solutions of equations of hydrodynamics. In the conformal case, up to 
second order in gradients, there are in total 5 available contributions consisting of the 3 hydrodynamic Weyl-invariant 
scalars multiplied by the velocity and 2 Weyl-invariant transverse vectors (|14l) 

J'' = i^^'"' + ^ + + (^i' ^1 + ^2 ^2 + i^3)«^} . (70) 

The overall factor of 1 /AGn in (j70p comes from the holographic expression for thermodynamic entropy. The on-shell 
divergencjfl of the current ()70|) was evaluated in reference [3| and reads, up to third order in gradients, 

cTpi, - hd{d - 2) [il - ^%^^^ ) w^AW^i. 



hd{d - 2) 



As understood in I38| for d — A, this expression makes it possible to constrain some of the coefficients appearing in 
([70|) . These arguments are based on the observation that local non-negativity should hold both when the shear tensor 
vanishes at a given point, as well as when it is arbitrary small (if it is large enough, then a^i^a^"^ dominates over other 
contributions and there are no further constraints). The first condition automatically implies that 

Jt^-^jl (72) 
whereas the second sets to zero all contributions which spoil non- negativity for very small cr^j/, i.e. 

jl+J^-lid-2)jl and .1=-^. (73) 



In the sense of conservation of the energy-momentum tensor given by II25I I 
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Note that appears in the divergence only in the combination j| + , so that shifting and jl keeping the sum 
constant does not change the divergence. This ambiguity comes from the freedom of modifying the entropy current by 
adding a multiple of the divergence- free term fe^-dp^^A^ = _^ g^^^^ Q, H 

and does not affect the 

local rate of entropy production. At third order there are no further constraints available so that j|' remains the only 
unspecified parameter afi^ecting the divergence ([7T|) . Results of [l^ 0| and 24 1 make it clear that j|' is not fixed by 
some higher order argument - dual gravitational constructions, which all guarantee non-negativity of the divergence, 
lead to different values of j|' . Thus, if the notion of local entropy production in the near-equilibrium regime makes 
sense, there must be some further constraints on the form of the hydrodynamic entropy current. This paper argues 
that one such constraint might be causality, which leads to considering the holographic entropy current based on the 
apparent horizon in the dual gravity description. 

The problem of constructing a candidate hydrodynamic entropy current on the gravity side of the correspondence 



251. These 



was first solved in 17] and then generalized to weakly curved boundary [30| and to arbitrary dimensions 
articles relied on using the bulk-boundary map defined by ingoing null geodesies supplemented with boundary diffeo- 
morphismj^l to map the area form of the black brane event horizon satisfying the area theorem onto a dual current 
of non-negative divergence. The main motivation for mapping bulk data along ingoing null geodesies was causality. 
Note however that such a constraint on the bulk-boundary map is self-consistent only when the bulk entropy carrier is 
causal^. Such a notion is provided by an apparent horizon, which along with the event horizon provides an example 
of a "generalized horizon" introduced in [24] . 

The geometric setup described in Section Hill contains a distinguished vector field v tangent to the horizon A. As 
anticipated in 2J] in the context of "generalized horizons" one motivation for introducing v is that the change of the 
area form on the horizon sections can be written in terms of the expansion 9 along v 

= -^^vVh, (74) 

where h is the determinant of the induced metric on the section. The generalized second law of thermodynamics is 
then the statement that the area of the leaves is non-decreasing under the above flow 



>0. (75) 



On the apparent horizon this area law is guaranteed by — 0, 0(„) < and C > 0. The boundary entropy current is 
obtained from v by means of rewriting the left hand side of (j75p within a chosen bulk-boundary map as a divergence 
of boundary current. This current is interpreted as a candidate boundary current and is given by [2^ 




where the prefactor involving Gn has been introduced to reproduce thermodynamic entropy at leading order and the 
AdS radius has been set to 1, as in ([T]). The technical assumptions used to derive (|76|) match those in Section Hill In 



The bulk-boundary map along ingoing null geodesies associates points on the apparent horizon, the event horizon or any other "gener- 
alized horizon" with boundary points lying on the same null geodesies moving close to the boundary in a direction specified by a given 
vector field. This vector field is taken to be proportional to in the leading order with subleading corrections modifying dual entropy 
current at orders higher than 2. As anticipated in section FlIII in the gauge JTTJ this bulk-boundary map acts trivially and maps points 
of the same x^^ position. Any such bulk-boundary map may be supplemented with boundary diffeomorphisms, which are generated 
by another vector field specified on the boundary. Such a vector field, if non-zero at leading order of the gradient expansion, must be 
also proportional to n^, which modifies the dual entropy current at second and higher orders. The only parameter in (I70I I shifted by 
boundary diffeomorphisms of such form is jf. For a detailed discussion of bulk-boundary maps see [itI . 

Relaxing the assumption of causality of bulk-boundary maps has so far not been explored. Note at this point that although the mapping 
along ingoing null geodesies seems (at least superficially) to be causal, boundary diffeomorphisms composed with a given bulk-boundary 
map might lead to causality violations (see Section |Vl for a discussion of this point). 
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particular, the formula (j76p is valid for a trivial bulk-boundary map, i.e. along null geodesies, which in the vicinity 
of the boundary move in the direction defined by u'^ . In the conformal case this direction can be modified only by 
second and higher order terms which change the entropy current at third order, and thus are beyond the scope of 
this article. The bulk-boundary map used here is not supplemented with boundary diffeomorphisms partly due to 
causality reasons (see Section |V] for more details). Because of this, the formula (1751) leads to the unique causal second 
order entropy current. 

To apply (|76p to the gravitational setup of Section one needs the form of v and a computation of the determinant 
of the bulk metric G. One finds, up to second order in gradients 

G = r2('*-i)5 (l - K,Si + ^7^2 • (77) 



V 2 (6r)2 , 

As discussed earlier, the vector v is completely fixed by the self-consistency of the bulk construction; the second order 
result (|54l) reads 

The right hand side of ([7S)) evaluated on the apparent horizon ([57]) leads to ([70)1 where ji and are fixed by 
Frobenius condition and equal 

^'1" " d{d-2)' = ^^^^ 
while the j|''s depend on the radial position of the apparent horizon and read 



d ' d{d-2)' 



,|l ^ (2 -3d) 
2d(d-2)' 

As a crosscheck one can easily see that the coefficients (l80t satisfy conditions (j72p and (|73)) . Comparing with the 
event horizon result 



2J, |25j one can see that the only difference is in the choice of j|' 



jl,EH ^ il,AH ^2 ■ ^^^^ 

Calculating the area theorem on the apparent horizon up to the second order in gradients one obtains 

^W-^^(")U.„ = ^> (82) 

which indeed matches the hydrodynamic result up to second order in gradients (I7ip . Calculating the third order 
contribution in the bulk requires third order geometry, which has so far not been obtained. Nevertheless the match 
is guaranteed by the formula ()76|) relating divergence of the entropy current to the area theorem on the apparent 
horizon modulo modifications of bulk-boundary map. 



V. SUMMARY 



This paper discusses the construction of apparent horizons in the geometry of conformal fluid-gravity duality. The 
motivation for this work are interrelated questions of local definition of entropy beyond equilibrium and foliation 
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dependence of apparent horizons of black holes. The reason for focusing on apparent horizons is that they are causal 
objects: they evolve in response to a flux of gravitational radiation or infall of matter. This makes them a preferable 
carrier for the notion of entropy beyond equilibrium in dual holographic field theories. There are however two caveats, 
which need to be taken into account: foliation dependence of apparent horizons [16| a nd locality of entropy production 
being directly related to the mapping of horizon information onto the boundary [l7[ • This paper focuses only on the 
first issue since the latter appears most severely at higher orders of the gradient expansion than are available for the 
geometry under consideration. 

The key idea behind this paper is that the apparent horizons of interest are only those which are covariant in the sense 
of the dual hydrodynamic description, i.e. can be covariantly specified (in the boundary sense) in terms of b, , g^i, 
and their gradients. This requirement chooses only those apparent horizons which have covariant dual hydrodynamic 
entropy currents. Apparent horizons which do not satisfy this condition evade a clear physical interpretation in terms 
of the dual field theory and are not studied in this paper. Because of the requirement of hydrodynamic covariance, 
this paper adopts the somewhat unusual strategy of first finding suitable null normals and only afterwards confirming 
that they are foliation- forming, rather than starting with a foliation and then proceeding to normals. The approach 
adopted here relies on the near-equilibrium regime, where one expects that at least one of the apparent horizons 
will "closely" follow the dynamics of the event horizon In the case of fluid-gravity duality this requirement is 
indeed satisfied - using even the results of 17[ alone one can easily check that the leading order event horizon is at 
the same time an (isolated) apparent horizon. It would certainly be interesting to try to apply similar methods to 
find apparent horizons in other black hole spacetimes, perhaps making contact with the framework of slowly evolving 
horizons 30-3^. 



The main result of this paper is that, up to second order in gradients in conformal fiuid-gravity duality, there exists a 
unique apparent horizon covariant in the hydrodynamic sense. It is very plausible that the uniqueness of this apparent 
horizon holds to all orders of the gradient expansion, as arguments in Section IIIII suggest. The apparent horizon in 
the geometry of fluid-gravity duality is isolated at leading and flrst subleading orders of the gradient expansion and 
becomes spatialonce second order gradient contributions are included. As expected 17| and confirmed by an explicit 
calculation in [24] , the apparent horizon gives rise to a notion of hydrodynamic entro py c urrent when the area form on 
the apparent horizon is mapped to the boundary in an appropriate way. Reference [17j introduced the map spanned 
along ingoing null geodesies, the main motivation for it being the causal structure of bulk spacetime. Such geodesies 
are specified by their tangent vector at the boundary and hydrodynamic covariance forces this to be proportional to 
the fluid velocity at leading order, but starting at second order additional contributions appear. These terms will 
modify the form of the entropy current at third and higher orders of the gradient expansion and are beyond the 
scope of this paper. The only freedom, which affects the divergence of an entropy current at second order comes from 
combining the bulk-boundary map specified by ingoing bulk geodesies with boundary diffeomorphisms 17]. It is not 
clear however, whether or not this leads to causality violations. One argument suggesting that it does was presented 



in 2J]. There, using the results from both Il7|,l25| and the present paper, it was shown that up to the second order 
in gradients the entropy current on the event horizon is equivalent to the entropy current on the apparent horizon 
when the bulk-boundary map in the latter case is supplemented with a particular boundary diffeomorphism. It would 
be very interesting to understand better the constraints on the form of the bulk-boundary map which follow from 
causality. 

It seems unlikely that explicit causality violations due to the choice of the bulk-boundary map can be visible at low 
orders of the gradient expansion in the same way as in the near-equilibrium regime it is hard to tell whether the 
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event or one of apparent horizons is a better entropy carrier [165 . Instead, one probably needs to look for some 
general principles or at concrete examples based on numerical solutions in asymptotically AdS spacetimes. One such 



a simple example is the gravity dual to boost- invariant flow [19|, |23|, |40 - 3] . In the boost-invariant case, the boundary 
dynamics depend on a single variable, the proper time r, and its large proper time limit is governed by boost-invariant 
hydrodynamics. Consider now the setup introduced in where the initial state given by the vacuum AdS space 
is excited in a boost-invariant way in the vicinity of r = by a time-dependent boundary metric. Such a quench 
leads to the emission of gravitational waves, and these propagate into the bulk and collapse forming a black hole. 
The black hole equilibrates and in the end is well-described in terms of fluid-gravity duality. In this example, the 
causal behavior of the apparent horizon is clearly visible only in the vacuum (before the quench) and in the far- 
from-equilibrium regime. However, close to equilibrium the apparent and event horizon follow each other closely as 



expected, and from that perspective there seems to be no reason to choose one over the other [16l |. Consider now a 
slight modification of the setup [l^ . As initial state at some late time ^ one can choose a boost-invariant black 
brane solution dual to perfect fluid hydrodynamics. Such a solution might have an arbitrary temperature and does not 
produce any entropy as required by perfect fluid hydrodynamics. Consider now the same kind of quench as considered 
in but now in the vicinity of r^. Such a quench will excite both far-from-equilibrium and hydrodynamics modes. 
The former equilibrate over a time scale set by the inverse of temperature, so by taking temperature to be large, one 
can effectively decouple them from analysis. Thus such setup serves as a causally clear example of entropy production, 
governed entirely by hydrodynamics. It would be very interesting to see what are the constraints on the part of the 
bulk-boundary map which corresponds to boundary diffeomorphisms following from this and similar examples. 

As for more obvious further projects, it would be very interestin g to calculate the location of apparent horizons in 
the cases of charged 47, 3|, non-conformal [49!], and superfluid 50, EH fluid-gravity dualities. In those examples 
there are more gradient terms available so that the relevant backgrounds may serve as further testing grounds for 
the claims of this paper. It would be also interesting to analyze the interplay between the position of the event and 



apparent horizons in the context of the AdS/CFT correspondence making use of the technology introduced in [21 [ 



In conclusion, the gravity dual to second order hydrodynamics of conformal media in arbitrary dimensions has a 
unique apparent horizon, which is covariant in the hydrodynamic sense. Possible ambiguities, which appear in the 
gravity construction, should not affect the amount of entropy produced between two equilibrium states (represented 
on the gravity side by two isolated horizons) . Furthermore if the foliation of the apparent horizon is fixed to all orders 
in the gradient expansion and there are very stringent constraints on the bulk-boundary map, then it is plausible that 
the local rate of entropy production in the near-equilibrium regime is a meaningful observable. This would be a very 
interesting result from the point of view of non-equilibrium statistical mechanics. 
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Appendix A: Weyl-covariance and available gradient terms 

The Weyl-covariant derivative V is defined so that for an arbitrary tensor Q'^ " of weight w (i.e. one that obeys 
Q'i " — >■ e'^'^Q^'" under Weyl transformations (O), Vx Q'^ " transforms homogeneously with weight w 

Vx Qi::, = Va Qi::, + w AxQ^,::. 

+ [.9Aa^^ - ^Ao, - <5^:^a] qZ::. + . . . (Ai) 
- - biA. - bZAx\ Q^::: - . . . 

The field A^l is given by ([5]) and transforms as a connection under Weyl transformations ([5]) . 

At first order in gradients there are two available contributions: 'DpUi, and 'Dph. The latter quantity vanishes at this 



order of the gradient expansion when evaluated on solutions of the equations of perfect fluid hydrodynamics |29j and 
thus at the first order the only nontrivial contribution comes from V^u,^. Taking its symmetric part one obtains the 
shear tensor a^^ 

1 



whereas its antisymmetric part is the vorticity of the flow 

1 



= A^^.]- (A3) 



Heve A(^^-) = + denotes symmetrization, whereas A[^^] = A^j^ — Aj^^ antisymmetrization. One can show 
that both tensors (jA2p and (|A3|) are transverse, i.e. u^^cr^u — it^a;^^ = 0. 



Following [25] one defines the Weyl-covariant Riemann tensor TZ/^i^xa 

'R-l^iyXcr = Ri^iyXcr + Vf^yl^^J^A^T - S[^9^][\'>'^] aAfj + AaAp - ^ffc/J^ (A4) 

Note that the Weyl-covariantized curvature tensors do not vanish even if the fluid lives in a fiat spacetime. With this, 
the Weyl-covariant Ricci tensor Tl^„ and Weyl-covariant Ricci scalar TZ appearing in (fT3)) can be defined: 



= n^xu^ and 7^ = 7^A^ (A5) 
Finally, the Weyl curvature tensor entering (jl5p is given by 

C^i.A<T = Rf.i'Xa + {tI^^u - 2{d-l) ) ■ ^^^^ 
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Appendix B: Example: the Vaidya metric 



To see how the method for finding apparent horizon introduced in the context of "generalized horizons" in 2J] and 
reviewed in Section [IIII works, consider the example of the Vaidya metric!^ 

ds^ = 2dwdr ~ f{r, w)dw^ + r^{de^ + sin^ Orcf^) (Bl) 

with 

^1-^. (B2) 
r 

In this case symmetry suggests taking S ~ r — g{w) (i.e. spherically symmetric apparent horizon). Then the covector 
m is 

rUa = {dr - g'{w)dw)^ (B3) 

and raising the index leads to 

m'^ ^{d^ + {f{r,w)-g'{w))drT . (B4) 

Here latin indices are used for {r,w,9,(p) and greek ones (below) for {w,d,(j)) - this is sort of analogous to the 
conventions used earlier for the AdS case. Choosing slicing on the horizon given by = {8^)^ leads to 

v- = {d^,+g'{w)drT (B5) 

Using gives to 

l = d.^ + ^f{r,w) dr and n ^ Xif{r,w) - 2g'{w)) dr. (B6) 
Imposing the normalization condition (1501) to fix A gives the expected result (s^ 

I — + 2/(^1 ^) and n = —dr- (B7) 
Note that in spherical symmetry vectors I and n are fixed up to an overall scaling by conditions 

W = Ubu'' = and Ian" = -1. (B8) 



The letter w denotes the ingoing Eddington-Finkelstein time coordinate to avoid clashing with the choice of v for the vector which 
defines the slicing 



